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We investigate the spectroscopic fingerprints of electron-plasmon coupling in integrated (PES) and 
angle-resolved photoemission spectroscopy (ARPES). To account for electron-plasmon interactions 
at a reduced computational cost, we derived the plasmonic polaron model, an approach based on the 
cumulant expansion of many-body perturbation theory that circumvents the calculation of the GW 
self-energy. Through the plasmonic polaron model, we predict the complete spectral functions and 
the effects of electron-plasmon coupling for Si, Ge, GaAs, and diamond. Si, Ge, and GaAs exhibit 
well defined plasmonic polaron band structures, i.e., broadened replica of the valence bands red- 
shifted by the plasmon energy. Based on these results, (i) we assign the structures of the plasmon 
satellite of silicon (as revealed by PES experiments) to plasmonic Van Hove singularities occurring 
at the L , D, and X high-symmetry points and (ii) we predict the ARPES signatures of electron- 
plasmon coupling for Si, Ge, GaAs, and diamond. Overall, the concept of plasmonic polaron bands 
emerges as a new paradigm for the interpretation of electronic processes in condensed matter, and 
the theoretical approach presented here provides a computationally affordable tool to explore its 
effect in a broad set of materials. 


I. INTRODUCTION 

The coupling of electrons and bosonic excitations, such 
as phonons, magnons, and collective charge-density fluc¬ 
tuations ( plasmons ), may significantly alter the elec¬ 
tronic properties of solids, triggering emergent phe¬ 
nomena which may not be explained within the con¬ 
text of a simple single-particle picture. A prototypi¬ 
cal example is the formation of Cooper pairs, 1 result¬ 
ing from electron-phonon coupling, and the emergence of 
superconductivity. 2 Whereas phonons have a characteris¬ 
tic energy in the range of 10—100 meV, the energy of plas¬ 
mons is typically between 1 and 10 eV. The spectral sig¬ 
natures of electron-plasmon coupling, thus, occur at en¬ 
ergies comparable to those of the ordinary quasiparticle 
states which are easily accessible in integrated (PES) and 
angle-resolved photoemission spectroscopy (ARPES). 

Theoretical approaches for the description of electron- 
plasmon interactions from first principles are generally 
based on many-body perturbation theory by combining 
the one-shot GW approximation 3,4 with the cumulant 
expansion approach 5,6 (GW+G). The cumulant expan¬ 
sion, formally exact for the case of an isolated core¬ 
electron, 6 stems from the solution of an electron-boson 
coupling Hamiltonian analogous to the case of electron- 
phonon coupling. The applicability of this method 
has been extended beyond the case of core satellites 
and its has proven useful in the description of quasi¬ 
particle dynamics, 7,8 and valence plasmon satellites in 
the integrated photoemission spectra of the homoge¬ 
neous electron gas, 9,10 simple metals, 11 silicon, 12 and 
graphene. 13,14 

First-principles calculations based on the 
GVF+cumulant approach have recently unveiled novel 
spectroscopic signatures of electron-plasmon coupling 
in the angle-resolved photoemission spectra (ARPES) 
of solids: the band structures of plasmonic polarons. 15 
Plasmonic polarons are elementary excitations resulting 


from the simultaneous excitation of an electron, and a 
plasmon and lead to the formation of dispersive spectral 
features in systems characterized by well defined plas¬ 
monic excitations (as, e.g., silicon, and transition-metal 
dichalcogenides) which follow closely the dispersive char¬ 
acter of the ordinary quasiparticle band structure. In 
particular, the dispersive nature of plasmonic polarons 
leads to the formation of well-defined plasmonic polaron 
bands which manifest themselves as broadened replica 
of the ordinary quasiparticle band structure, red-shifted 
by the plasmon energy. The existence of plasmonic 
polaron bands has recently been confirmed by ARPES 
measurements of the full band structure of silicon in 
an energy window that extends up to 35 eV below the 
Fermi energy. 16 

In this work, we introduce the plasmonic polaron 
model , an approach for the computation of plasmonic po¬ 
laron band structures that circumvents the calculations 
of the GW self-energy. Beside illustrating analytically 
the emergence of plasmonic polaron bands within the 
GW+cumulant formalism, the plasmonic polaron mod¬ 
els allows to account for the effects of electron-plasmon 
coupling in the spectral properties of solids at a reduced 
numerical cost. We thus discuss electron-plasmon inter¬ 
actions in Si, GaAs, Ge, and diamond. Our results reveal 
that the different substructures of the plasmon satellite 
of silicon, as measured by XPS, may be assigned to Van 
Hove singularities arising from the plasmonic polaron 
bands taking place at the L, O, and X high-symmetry 
points of the first Brillouin zone. For GaAs, Ge, and 
diamond, we predict plasmonic polaron bands in quali¬ 
tative agreement with silicon, and we identify the spec¬ 
tral features that may guide the observation of plasmonic 
polaron bands in these compounds. Overall our results 
indicate that plasmonic polaron bands provide a novel 
concept that may prove generally applicable in the inter¬ 
pretation of the integrated and angle-resolved spectral 
properties of solids. 
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II. PLASMONIC POLARON MODEL 

In the following, we derive an approach to predict 
the band structures of plasmonic polarons at the nu¬ 
merical cost of a band structure calculation, circum¬ 
venting the explicit calculation of the GW self-energy. 
The cumulant expansion for the single-particle Green 
function 6,9-14,17,18 provides an ideal framework for the 
description of electron-plasmon interactions in solids. In 
particular, the starting point for the following discussion 
is provided by the derivation of the GW +cumulant ap¬ 
proach reported in Ref. 11 (referred to as GW + Cahk 
hereafter). If one neglects rare events in which multiple 
plasmons are emitted simultaneously with the excitation 
of a photo-hole, the GW+Gark spectral function can be 
expressed as : 11 

A(k, w ) = 5>Q p (k, w ) + AQ p (k. w ) * k jW )]. (1) 

n 

The first term in Eq. (1) accounts for the spectral signa¬ 
tures of quasiparticle excitations (i.e., in absence plasmon 
emission) and is defined as: 

_r rak (o;)_ 

7 r [w - e nk - AE nk (e nk )] 2 + [r„ k (w )] 2 ’ 

( 2 ) 

where AE„ k denotes the GoWo quasiparticle 
correction 3,4 to the Kohn-Sham eigenvalues 19,20 e nk , and 
r nk = l/r„ k the quasiparticle lifetime. Here and in the 
following we adopt Hartree atomic units. As discussed 
in Ref. 15, the convolution product in the second term 
of Eq. (1) introduces novel dispersive features in the 
spectral function, which account for events in which a 
photo-hole and a plasmon are excited simultaneously. 
The term is defined as : 11 
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where /3 nk (uj) = 7 r _ 1 ImE nk (e nk - u)9(^ - u). 

To simplify the evaluation of Eq. (1), we notice that the 
self-energy of systems characterized by well-defined plas¬ 
mon resonances - as, e.g., simple metals 11 (A1 and Na) 
and semiconductors 21 (Si and GaAs) - exhibits an imag¬ 
inary part with a large plasmon peak at w = £ rak — 
where fl p is the plasmon energy. The plasmon peak 
in ImE stems directly from the plasmonic resonance of 
Ime _ 1 (u;), with e being the dielectric function in the 
random-phase approximation (RPA). Therefore, the sec¬ 
ond term of Eq. (1) may be simplified by assuming 
ImE„ k ~ A nk L(w), where A nk is a constant and L a nor¬ 
malized Lorentzian function centered at u = £ rak — fl p , 
that is: L(w) = (j/tt )/[(w - £ nk + H p ) 2 + q 2 ]. The 
ansatz introduced above for ImE is an oversimplification 
of the structures of the self-energy that, as demonstrated 
below, suffices to reproduce the essential experimental 
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Figure 1. Angle-resolved spectral function of Silicon for mo¬ 
menta along the I -A direction evaluated from (a) the plas¬ 
monic polaron model [Eq. (5)] and (b) converged SGVF + 
Cahk calculations from Ref. 15. The vertical arrow indicates 
the plasmon energy Q p of silicon. The bare PBE band struc¬ 
ture and its shifted plasmonic polaron replica are shown in 
blue. 



spectral features arising from the electron-plasmon inter¬ 
action. 

The Lorentzian ansatz for ImE can be derived 
within the framework of the ordinary plasmon-pole 
approximation 4,22-28 for GW calculations by considering 
non-dispersive bands and non-dispersive plasmons. 

If the full-width at half maximum (FWHM) 7 is as¬ 
sumed to be sufficiently small (7 -C |£ nk — H p |) one may 
retain only the first term in the numerator of Eq. (3), 
which simplifies to: 
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The convolution product in Eq. (1) can be approximated 
using Eq. (4) and the spectral function reduces to: 


A(k,w) ~ 

n 
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Equation (5) - the central equation of the plasmonic po¬ 
laron model illustrates the emergence of band struc¬ 
ture replicas as a result of electron-plasmon interactions. 
The resulting additional spectral features are the band 
structures of plasmonic polarons - copies of the ordinary 
quasiparticle bands red-shifted by the plasmon energy 
f l p . These bands, which follow closely the momentum 
dependence of the quasiparticle energies, reveal the dis¬ 
persive nature of plasmonic polarons. 

In Eq. (5), the prefactor A, lk / 7 rfl 2 determines the rel¬ 
ative spectral weight of the plasmonic polaron features 
with respect to the quasiparticle peaks, and it may be 
related to the quasiparticle weight ,Z„ k . In practice, the 
quasiparticle part of the spectral function is related to 
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Figure 2. (a) Spectral function of silicon determined from the plasmonic polaron model [Eq. (5)]. The plasmonic polaron 
intensity has been magnified by a factor of 3 to enhance visibility, (b) Experimental (integrated) x-ray photoemission spectrum 
(XPS) of silicon reproduced from Ref. 18. (c) Quasiparticle and plasmonic polaron (magnified by a factor of 2) density of states 
(DOS). Different peaks and shoulders of the plasmonic polaron resonance in the XPS spectrum are attributed to plasmonic 
Van Hove singularities occurring on different high-symmetry lines in the Brillouin zone (red circles). 


the quasiparticle weight through (aj)doj = Z n k- 

Similarly, making use of Eq. (4) and of the normalization 
of L(w) one can write: 

j [^ P (k,w) *A%(k,u})\duj = ^Z n k. (6) 

Requiring the normalization of the spectral function 
(that is, A n k(u) = 1) it is possible to express A n k 
explicitly in terms of Q p and Z n k' 29 

A -k= (7) 

To further simplify the calculation of the plasmonic po¬ 
laron band structures, we neglect the quasiparticle cor¬ 
rection to the Kohn-Sham eigenvalues (i.e., AE„k = 0) 
in Eq. (2), and we introduce a simple analytical model 
for the quasiparticle lifetime. According to Fermi liquid 
theory, the inverse lifetime of quasiparticle states close 
to the Fermi energy increases quadratically with their 
energy , that is r~ k oc (e n k — p) 2 ■ Here we consider 
semiconductors, whereby ImE(w) = 0 for energies in the 
range w £ [fi — E g , p ], E g being the quasiparticle band 
gap, owing to the absence of (electronic) decay channels 
that may induce a de-excitation of the photohole. We 
combine these two ideas to define a simplified model for 
the electronic lifetimes of semiconductors. Henceforth, 
we assume a quadratic dependence of the quasiparticle 
lifetimes on the frequency (relative to the Fermi energy), 
in the form: T(w) = p + a(w + E g ) 2 6(p — E g — w). In 
short, rj = 0.05 eV is introduced to avoid divergences 


close to the Fermi energy, whereas the term E g ensures 
that the inverse lifetime of states in the energy window 
lj £ [p — E g , fj] is negligible due to the absence of decay 
channels (phonon-assisted decay processes are neglected 
here). The parameter a is determined by imposing the 
equality between the integrated area of the broadening 
model defined above and the area underlying the imagi¬ 
nary part of the self-energy, which coincides with A„k- In 
other words, we require jff T(w)dw = jff ImE„k(w)dw = 
A n k, where w denotes the energy such that ImE„k(w) cs 0 
for w < w [in the following, w = 2(e„k — f2 p )]. Alterna¬ 
tive broadening function models for T may be employed, 
however the precise choice for T does not alter the quali¬ 
tative features of the plasmonic polaron band structures 
determined through Eq. (5). 

In summary, Eqs. (2) and (5) define a simple gen¬ 
eral procedure for predicting the plasmonic polaron band 
structures of systems characterized by well defined plas¬ 
monic excitations, avoiding the explicit calculation of the 
GW self-energy. The resulting model requires knowledge 
of a few parameters which, for semiconductors, may eas¬ 
ily be inferred from experiment and/or first-principles 
calculations: the plasmon energy Q p , the quasiparticle 
weight Z n k, and the quasiparticle band gap E g . For sp- 
bonded semiconductors typically 0.7 < Z n k < 0.9. In the 
following, we assume Z n k — 0.8, as, e.g., for bulk Si and 
Ge. To further reduce the number of external param¬ 
eters, we obtain the plasmon energy Q p from a simple 
model for the frequency-dependent dielectric function of 
semiconductors: 30 Q 2 = w 2 + E 2 , where w p is the ho¬ 
mogeneous electron gas plasma frequency at the valence 
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Figure 3. Spectral function of (a) Ge and (c) GaAs determined from the plasmonic polaron model [Eq. (5)]. The plasmonic 
polaron intensity has been magnified by a factor of 3 to enhance visibility, a, /?, and 7 denote regions characterized by low 
density of states due to the absence of plasmonic polarons and quasiparticle states. Panel (b) and (d): quasiparticle and 
plasmonic polaron DOS of Ge and GaAs. 
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electron density 30 and E g the quasiparticle band gap. In 
the following, we consider the experimental value for the 
quasiparticle band gap E g . Once the parameters E g , 
and Z are given, the evaluation of Eq. (5) can be per¬ 
formed through a simple band-structure calculation. In 
the following, all calculations are performed in a plane- 
wave basis within the Quantum Espresso code. 31 ' 32 

Before moving on to discuss electron-plasmon interac¬ 
tions in semiconductors, we briefly summarize the ap¬ 
proximations employed in the plasmonic polaron model 
and we comment on their validity. The only assumptions 
used until this point are that (i) the imaginary part of 
the self-energy can be approximated by a sharp pole at 
the plasmon energy and (ii) the photo-emission process 
can be described within the sudden approximation. 

The assumption (i) is based on the fact that systems 
characterized by plasmonic excitations typically exhibit 
a well defined resonance in the imaginary part of the 
inverse dielectric function e (owing to its relation to 
the electron energy loss function 33 ) which introduces a 
corresponding plasmonic resonance in ImE via the re¬ 
lation ImE = Im[e -1 t>]G. This approximation is thus 
justified for systems in which the plasmon-pole approx¬ 
imation yields a reasonable description of the inverse 
dielectric matrix. The plasmon energy dispersion, ne¬ 
glected within approximation (i), would introduce an ad¬ 
ditional broadening of the self-energy but it is not ex¬ 
pected to alter the qualitative spectral features obtained 
from the model. The approximation (ii) in practice cor¬ 
responds to neglecting extrinsic effects to the photoemis¬ 
sion process 12,34-36 . Extrinsic effects, which account for 
the interactions between the photoelectron and the sys¬ 
tem after emission, are expected to introduce an addi¬ 
tional broadening of the plasmonic polaron features and 
a renormalization of their intensity 12 ’ 18 . Additionally, in 
the following we estimate the plasmon energies from a 


simplified dielectric function model for semiconductors 
which assumes that collective charge density fluctuations 
involve the entirety of the valence electrons and that elec¬ 
tronic states are sufficiently delocalized 37,38 (conditions 
obeyed by sp-bonded systems as those considered here). 
For d- and /-electron systems, where these conditions 
may not apply, it would be more appropriate to employ 
plasmon energies derived either from experiment, from 
first-principles calculations of the energy-loss function in 
the random-phase approximation, or from more elabo¬ 
rate dielectric function models. 


III. PLASMONIC POLARONS AND 
PLASMONIC VAN HOVE SINGULARITIES 

As a first step, we validate the plasmonic polaron 
model by comparing the spectral function of silicon 
obtained from Eq. (5) with accurate first-principles 
calculations based on the SternheimerGVF-|-cumulant 
approach 11 ’ 15,21 ’ 39 (SGIV + Gahk)- In short, SGIV pro¬ 
vides an accurate reference method for the calculation 
of the full spectral function. In particular, in SGW (i) 
summations over empty states are avoided by computing 
the screened Coulomb interaction W and Green’s func¬ 
tion through the iterative solution of the Sternheimer 
equation and (ii) beside the random-phase approxima¬ 
tion (RPA) there are no further approximations (as, e.g., 
the plasmon-pole model) involved in the computation of 
the dielectric function. More details on the theoretical 
and numerical aspects underlying the Sternheimer-GIV 
approach can be found elsewhere. 21,39 ’ 40 In Fig. 1, we 
report the angle-resolved spectral function of Silicon for 
momenta within the first Brillouin zone along the T-A 
high-symmetry line as obtained from (a) the plasmonic 
polaron model [Eq. (5)] and (b) the SGIV + Gahk ap- 
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proach. The PBE band structure of silicon and its plas- 
monic polaron replica (red-shifted by the plasmon energy 
f l p ) are superimposed as a continuous blue line for com¬ 
parison. Overall, as illustrated in Fig. 1, the proposed 
plasmonic polaron model allows us to reproduce the qual¬ 
itative features of SGW + Cahk approach at the cost of 
a band-structure calculation. Based on this finding, we 
now move on to discuss the full spectral function of Si, 
GaAs, Ge, and diamond and the spectral fingerprints of 
plasmonic polaron excitations in these compounds. 

In Fig. 2 we report (a) the full spectral function of sil¬ 
icon evaluated within the plasmonic polaron model and 
(c) the corresponding density of states (DOS). The plas¬ 
monic polaron DOS has been calculated by ignoring the 
frequency dependent component of the broadening func¬ 
tion (a = 0) to emphasize the structure arising from 
the different Van Hove singularities. To approximately 
account for the different cross-section effects, the DOS 
has been obtained as a weighted sum of the s-orbital 
and p-orbital projected-DOS, with weights given by the 
relative cross-section of the s and p states at a phonon 
energy of 800 eV 41 . The plasmonic polaron DOS has 
been magnified by a factor of two to enhance the visi¬ 
bility on the same scale of the quasiparticle spectral fea¬ 
tures. The experimental integrated x-ray photoemission 
spectrum (XPS) reproduced from Ref. 18 is reported for 
comparison in Fig. 2 (b). The spectral function exhibits 
a set of bright quasiparticle bands for binding energies 
in the range of 0 to 12 eV. These bands arise from the 
first term of Eq. (5) and define the ordinary band struc¬ 
ture of silicon. Their intensity is proportional to the 
emission rate of a photoelectron in absence of plasmonic 
excitations. In addition to the ordinary quasiparticle 
band structure, electron-plasmon coupling introduces ad¬ 
ditional dispersive spectral features in the angle-resolved 
spectrum. These spectral features, which manifest them¬ 
selves as broadened replica of the valence band structure, 
reveal the excitation of plasmonic polarons, i.e., elemen¬ 
tary excitation in which the energy of the absorbed pho¬ 
ton contributes to the emission of a photoelectron and 
the excitation of a plasmon with energy ~ Q p . 

Our previous work 15 revealed that, similarly to ordi¬ 
nary quasiparticle bands, also plasmonic polaron band 
structures may lead to the formation of Van Hove singu¬ 
larities in the density of states (DOS). For quasiparticle 
bands, the density of states J may be expressed as: 

<«> 

where 5k denotes the surface e n k = constant. Equa¬ 
tion (8) indicates that whenever the first momentum 
derivative of the quasiparticle bands e n k vanishes (that 
is, |Vk£nk| = 0) the density of states exhibits sharp reso¬ 
nances, known as Van Hove singularities. 42 Consequently, 
the structures that characterize the DOS for binding en¬ 
ergies between 0 and 12 eV may be attributed to specific 
high-symmetry points of the first BZ, whereby quasipar¬ 
ticle bands are flat, and thus their derivative vanishes. 
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Figure 4. Spectral function (a) and DOS (b) of diamond 
determined from the plasmonic polaron model [Eq. (5)]. 

The concept of Van Hove singularities may easily be 
generalized to the case of plasmonic polarons. In this 
case, however, we do not expect the emergence of sharp 
resonances owing to the large broadening characteristic 
of plasmonic polaron excitations. For silicon the XPS 
measurements at binding energies between 15 to 30 eV 
reveal a broad plasmon satellite characterized by a sub¬ 
structure of peaks and shoulders. These structure can 
be assigned to Van Hove singularities that arise from the 
plasmonic polaron band structures. 

Having established a numerically efficient procedure 
to compute the plasmonic polaron contribution to the 
full spectral function, we can now proceed to assign the 
structures in the XPS plasmon satellite of silicon to high- 
symmetry points in the first BZ. As illustrated in Fig. 2, 
the two main peaks that characterize the plasmon satel¬ 
lite at 20 and 24 eV are Van Hove singularities result¬ 
ing from a flattening of the plasmonic polaron bands at 
X and H (the middle point on the W-L high-symmetry 
line), respectively. Similarly, the two shoulders at 18 and 
26 eV originate from the L point. The DOS of plas¬ 
monic polarons reported in Fig. 2 (c) further validates 
this analysis, illustrating that the Van Hove singularities 
of the plasmonic polaron bands leads to spectral features 
in excellent agreement with the experimental reference 
data. These results provide a simple rationale to unravel 
the complex features of the plasmon satellites in the in¬ 
tegrated photoemission measurements of solids. In par¬ 
ticular, plasmon satellites can be attributed to plasmonic 
Van Hove singularities and their substructures arise from 
the flattening of the plasmonic polaron bands at different 
high-symmetry points in the Brillouin zone. 

We now move on to discuss the plasmonic polaron 
band structures of Ge and GaAs. In Fig. 3 we report 
the full spectral function of Ge (a) and GaAs (c) deter¬ 
mined from the plasmonic polaron model and the corre¬ 
sponding DOS in panels (b) and (d), respectively. The 
spectral functions of Ge and GaAs exhibit similar qual¬ 
itative features to Si. For binding energies up to 15 eV 
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below the Fermi energy, the spectrum is characterized 
by the ordinary quasiparticlc band structure. For bind¬ 
ing energies larger than 15 eV, the spectral function is 
dominated by the plasmonic polaron bands. These bands 
are expected to provide a dominant contribution to the 
density of states in this energy range. Overall for Si, 
Ge, and GaAs, the broadening of the plasmonic polaron 
bands may obstruct the experimental observation of the 
individual bands. However, based on the results pre¬ 
sented in Figs. 2 and 3, we identify general patterns that 
may drive the experimental observation of plasmonic po¬ 
laron bands. In particular, Si, Ge, and GaAs manifest 
a depletion of density of states around the T point for 
binding energies between 20 and 30 eV. As illustrated in 
Fig. 3, this leads to the emergence of a diamond-shaped 
structure (denoted as 7 in Fig. 3) in ARPES measure¬ 
ments. The low-intensity region 7 can be attributed to 
the gap between the first three degenerate plasmonic po¬ 
laron bands (that is, the plasmonic replica of the three 
top-most quasiparticle bands) at the T point and the 
fourth plasmonic polaron band. Similarly, we observe 
a lowering of spectral intensity around W and X (de¬ 
noted as a and ft in Fig. 3) for binding energies between 
11 and 18 eV below the Fermi energy. In this case, the 
a and ft regions stem from the gap between the lowest 
quasiparticle band and the first plasmonic polaron band. 

Our calculations for the plasmonic polarons DOS 
[Fig. 3 (b) and (d)] indicates that, similarly to the quasi¬ 
particle DOS, plasmonic polarons leads to the formation 
of well defined Van Hove singularities. We thus predict 
that XPS measurements of Ge (GaAs) should reveal a 
plasmon satellite characterized by at least three different 
structures at binding energies of 18.5, 23.0, and 26.1 eV 
(18.1, 22.0, and 26.3 eV). As in the case of silicon, we 
attributed these structures to the flattening of the plas¬ 
monic polaron bands at different high-symmetry points 
of the Brillouin zone (the X , f2, and L points) as illus¬ 
trated in panels (b) and (d) of Fig. 3. 

As a last example, we consider the case of diamond. 
The electronic properties of diamond are qualitatively 
different from those of the semiconductors discussed 
above. In particular, diamond is wide gap insulator 
(Eg ~ 5.5 eV), its bands are more dispersive than Si, 
Ge, and GaAs, and it has a larger plasmon energy 
il p ~ 22.5 eV. These features are expected to affect con¬ 
siderably the band structures of plasmonic polarons. An 
inspection of the full spectral function [Fig. 4 (a)] reveals 
that the plasmonic polaron bands emerge at larger bind¬ 
ing energies as compared to the other compounds, owing 
to the larger plasmon energy Because of the shorter 
lifetime of plasmonic polarons at these binding energies, 
however, their spectral signatures appear broader and 
less intense. For the case of diamond, electron-plasmon 
interactions are thus expected to introduce low inten¬ 
sity spectral features that might be difficult to reveal in 
ARPES measurements. However, a possible evidence of 
electron-plasmon coupling might arise in the form of an 
V-shaped resonance at the bottom of the valence band 


[denoted as 7 in Fig. 4 (a)] which results from the merging 
of the lowest-energy valence band and the highest-energy 
plasmonic polaron band at the T point. Similarly, the 
large bandwidth of quasiparticle and plasmonic polaron 
bands of diamonds leads to less pronounced Van Hove 
singularities as compared to silicon and GaAs [Fig. 4 (b)]. 
Correspondingly, the plasmon satellite is unlikely to ex¬ 
hibit substructures that may be attributed to Van Hove 
singularities occurring at different points in the Brillouin 
zone. In particular, the broadening is expected to wash 
out and merge the two resonances at 34.2 and 38.2 eV 
indicated in Fig. 4 (b). 


IV. CONCLUSIONS 

In summary, we presented a study of the spectroscopic 
signatures of clectronic-plasmon coupling in integrated 
and angle-resolved photoemission experiments on tetra¬ 
hedral semiconductors. Based on the cumulant expan¬ 
sion for the spectral function, we derived the plasmonic 
polaron model, a simplified approach for the calculation 
of the band structures of plasmonic polarons in solids. 
This model, that does not require the explicit compu¬ 
tation of the GW self-energy, allows us to predict the 
full spectral function of solids at the numerical cost of 
a band structure calculation. Our results unveil the ef¬ 
fects of electron-plasmon coupling for Si, Ge, GaAs, and 
diamond. Si, Ge, and GaAs exhibit reasonably well de¬ 
fined plasmonic polaron band structures, that manifest 
themselves as broadened replica of the valence bands red- 
shifted by the plasmon energy, whereas for diamond plas¬ 
monic polarons spectral features are suppressed by short 
lifetime effects. Additionally, this study reveals that the 
different structures that characterize the plasmon satel¬ 
lite of silicon may be attributed to plasmonic Van Hove 
singularities occurring at specific high-symmetry points 
in the Brillouin zone. This indicates that the dispersive 
nature of plasmonic polarons may be revealed already in 
integrated photoemission experiments. Finally, our anal¬ 
ysis provides clear guidelines for identifying the spectral 
signatures of plasmonic polarons in ARPES. In partic¬ 
ular, for Si, Ge, and GaAs, plasmonic polaron bands 
should lead to low-intensity regions in the spectral func¬ 
tion around the X and T points for binding energies 
around 15 and 25 eV, respectively. 

Overall, plasmonic polaron bands emerge as an impor¬ 
tant new concept for the interpretation of integrated and 
angle-resolved photoemission experiments. These find¬ 
ings calls for renewed experimental efforts along the lines 
of Ref. 16 for silicon. The work presented here provides 
a simple predictive tool that may prove useful to unravel 
the complexity of photoemission spectra for a broad class 
of materials (including, e.g., semiconductors, metals, and 
^//-electron systems), and to advance our understanding 
of elementary electronic excitations in condensed matter. 
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